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Abstract 

We construct an event-based computer simulation model of the Einstein-Podolsky-Rosen-Bohm 
experiments with photons. The algorithm is a one-to-one copy of the data gathering and analysis 
procedures used in real laboratory experiments. We consider two types of experiments, those with 
a source emitting photons with opposite but otherwise unpredictable polarization and those with 
a source emitting photons with fixed polarization. In the simulation, the choice of the direction of 
polarization measurement for each detection event is arbitrary. We use three different procedures 
to identify pairs of photons and compute the frequency of coincidences by analyzing experimental 
data and simulation data. The model strictly satisfies Einstein's criteria of local causality, does not 
rely on any concept of quantum theory and reproduces the results of quantum theory for both types 
of experiments. We give a rigorous proof that the probabilistic description of the simulation model 
yields the quantum theoretical expressions for the single- and two-particle expectation values. 
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I. INTRODUCTION 



A fundamental problem, originating from the work of Einstein, Podolsky and Rosen 
(EPR) [1] and reformulated by Bohm ^ is to explain how individual detection events, 
registered by different detectors in such a way that a measurement on one particle does 
not have a causal effect on the result of the measurement on another particle (Einstein's 
criterion of local causality), give rise to the two-particle quantum correlations that are found 
in experiments 3], 4, 5, 6, 3, S, S H, 11, 12]. 

In Einstein-Podolsky-Rosen-Bohm (EPRB) experiments, individual events are registered, 
correlations between them are calculated and are found to correspond to the two-particle 
correlation for the singlet state. Since, in quantum theory, the basic equation that describes 



individual events is not known 



131 ] . quantum theory simply cannot be used to construct 



a numerical algorithm to simulate the individual events. Of course, using pseudo-random 
numbers we could generate events according to the probability distribution that is obtained 
by solving the Schrodinger equation. However, the challenge is to explain how the individual 
events can give rise to the two-particle correlations of the singlet state without invoking 
concepts of quantum theory. 

The question that we address in this paper is: Given the existing experimental data 
(numbers recorded during an experiment, stored on computer disks, and analyzed long after 
the data is taken), that, when analyzed properly, yields expectation values which are in 



good agreement with the predictions of quantum theory a ay, ay , |l2J], is it possible 
to construct an event-based simulation algorithm that satisfies Einstein's criteria of local 
causality, generates the same kind of data as in experiment, and is capable of reproducing 
exactly the single- and two-particle averages of quantum theory for a system of two S — 1/2 
particles? 

Within the context of local realist probabilistic models, a rigorous proof that the existence 
of an algorithm that describes the outcome of real EPRB experiments cannot be excluded, 
has been given earlier |14]. Although such a proof is very valuable, actually finding such 
algorithms using local, causal processes to generate the probability distributions of quantum 
theory, is another challenge. 

In this paper, we present results of a complete simulation of Aspect-type experiments 
using an Einstein local, causal event-based simulation model. An important feature in these 
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experiments is the arbitrariness in the choice of the directions in which the polarization will 
be measured, for each individual detection event , 0]. This feature has not 

been taken into account in our earlier work [ijj] but is fully accounted for in the simulation 
procedure that we describe in this paper. 

The paper is organized as follows. In Section [III we describe the experimental set-up, 
the data gathering method and the data analysis procedures used in EPRB experiments 
with photons, closely following Ref. jl], The sources used in EPRB experiments emit 
photons with opposite but otherwise unpredictable polarization. Each photon propagates 
to an observation station consisting of a polarizer and two detectors. In accordance with 
quantum theory and experiments, we expect the two-particle correlation to agree with the 
expression obtained by assuming that the quantum state is a singlet. We refer to this 
experimental set-up as Case I. Inserting polarizers between the source and the observation 
stations changes the pair generation procedure such that the two photons that enter the 
observation stations have a fixed polarization. In this case, the photon intensity recorded by 
the detectors behind the polarizers in each observation station obeys Malus law. We refer 
to this set-up as Case II. A brief review of the quantum theoretical description of Cases I 
and II is given in Section IIII1 

In real experiments, macroscopic or microscopic, we need a well-defined procedure to 
decide if two detection events stem from a single system. In real EPRB experiments with 
photons, the time at which the events are registered is used for this purpose. However, the 
criterion that is used to select the events that stem from a single two-particle system is, to 
considerable extent, arbitrary. In Section IIVI we study this aspect by analyzing publicly 



available experimental data for an EPRB experiment with photons [8]. We present results 
of an analysis using three different procedures: 



First, we simply divide the time interval of measurement in equally spaced bins [16 ]. 
For each station, we determine the number of events per bin. From this data, we 
compute the coincidences. Effectively, this procedure compares the detection times at 
both stations with the time of a reference clock, using a coincidence window with a 
width that is equal to the bin size. 

Second, we employ the criterion used in the experiment We compute the 

coincidences of a detection event at station 1 and a detection event at station 2 by 
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comparing the time difference of these events with a fixed time window, that is we use 
relative times to determine the coincidences. 

• Finally, in the third procedure we first maximize the number of coincidences by 
shifting by the same amount, the detection times of station 2 relative to those of 
station 1, and then use the second procedure to count the coincidences. This two-step 
procedure reproduces the published results [8|. 

Our analysis shows that the first and third procedure may yield a result that is in reason- 
able agreement with the prediction of quantum theory if the bin size or coincidence window 
is sufficiently small. The data obtained by the second procedure is similar except that for a 



particular choice of the time window, the result is in conflict with quantum theory. 
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eral, these results support the idea that the idealized EPRB gedanken experiment 
that agrees with quantum theory cannot be performed 201 ] . 

In Section|V]we describe an Einstein local, causal event-based computer simulation model, 
based on the EPRB experiment with photons performed by Weihs et al . The crucial 



point of the present and of our earlier work [15 
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23] is that we simulate a model of 



the real EPRB experiments, not of the simplified, gedanken-type version that is commonly 
used [l8, Q. We give an explicit description of the algorithm to simulate the photons 
one by one, the observation stations containing the polarizers and detectors, and the data 
analysis procedure. The polarizers are modeled such that we reproduce the quantum theo- 
retical results for Case I and Case II without changing the algorithm for the polarizers, that 
is the functionality of all polarizers is the same. In contrast to the real EPRB experiments 
with photons [s, Q], the number of orientations per polarizer to choose from is not limited 
to two. 

Section IVII gives a rigorous analytical treatment of the probabilistic model of our simu- 
lation algorithm and proves that this model can reproduce the single-particle averages and 
the two-particle correlation of a system of two quantum spins for Cases I and II. This prob- 
abilistic model is identical to the one studied in Ref. |14]. except for the concrete model of 
the time-delay mechanism. In Section fVIII we present our simulation data and demonstrate 
that there is excellent agreement with the results obtained from quantum theory and the 
probabilistic model. In Section IVHH we study the effect of the time window on the fre- 
quency of coincidences and show that the simulation model readily reproduces published 
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FIG. 1: (color online) Case I: Schematic diagram of an EPRB experiment with randomly polarized 
particles. 

experimental data, including the statististics of the single-detection events. A summary of 
our results is given in Section HXl 

II. EPRB EXPERIMENT WITH PHOTONS 

A schematic diagram of Case I is shown in Fig. [H A source emits pairs of photons with 
opposite but otherwise unpredictable polarization. Each photon of a pair propagates to an 
observation station in which it is manipulated and detected. The two stations are separated 
spatially and temporally. This arrangement prevents the observation at station 1 (2) to have 
a causal effect on the data registered at station 2 (1). In Case II (see Fig. [2]), additional 
polarizers are inserted between the source and the observation stations [5| such that the 
two photons that enter the observation stations have a fixed polarization. We denote the 
orientations of these polarizers by the angles r/i and r)2- 

As the photon arrives at station z = 1,2, it passes through a polarizer. The orientation of 
the polarizer in observation station i is characterized by the angle 9,^ which may be chosen at 
random. As the photon leaves the polarizer, it generates a signal in one of the two detectors. 
Each station has its own clock (not shown) that assigns a time-tag to each signal generated 
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FIG. 2: (color online) Case II: Schematic diagram of an EPRB experiment with particles with 
fixed polarization. 

by one of the two detectors . Effectively, this procedure discretizes time in intervals, 
the width of which is determined by the time-tag resolution r. In experiment, the time-tag 
generators are synchronized before each run 8l |9|. This procedure is necessary because in 
time, the clocks may become unsynchronized 8|, [9(. 

Note that the description given earlier is only a pictorial description of real EPRB ex- 
periments with photons, as they are carried out in a laboratory. The experimental facts 
are the settings of the various apparatuses and the detection events. What happens in be- 
tween activating the source and the registration of the detection events is not, or cannot be, 
measured and is therefore not known. In this sense, the photon should be regarded as an 
element of a model or theory for the real laboratory experiment only. 

In the experiment, the firing of a detector is regarded as an event. At the nth event, 
the data recorded on a hard disk (not shown) at station i = 1, 2 consists of n> j, x n> i = ±1, 
specifying which of the two detectors behind the selected polarizer fired and the time tag 
t n) i indicating the time at which a detector fired. Hence, the set of data collected at station 
% — 1, 2 during a run of N events may be written as 

Tj {%n,ii tn,ii @n,i\Tl !>---> -A} • (1) 
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Any real EPRB experiment requires some criterion to decide which detection events are 
to be considered as stemming from a single two-particle system. In EPRB-type experiments 
with photons, this decision is taken on the basis of coincidence in time |8|, |24j]. However, 
as discussed in Section HVl this identification procedure is not unique. Coincidences can, 
for example, be identified by comparing the time differences {i nj i — t n ^\n = 1, . . . , N} with 



a window W 
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In this case, for each pair of rotation angles a and /?, the number of 
coincidences between detectors D x>1 (x = ±1) at station 1 and detectors D y 2 (y = ±1) at 
station 2 is given by 



C xy — C xy (a, 0) — J2n=l ^n,A,x n Ja,e n Jp,i 



■11.2 



xQ(W-\t n!l -t n!2 \), (2) 

where 0(t) is the Heaviside step function. The single-particle averages and correlation 
between the coincidence counts are then given by 



Ylx,y=±l X C X y 



Ylx,y=±l 1)Cxy 



Ei,j/=±1 ^xy 



Ei, s =±i x yC xy 



^2x,y=±l Cxy 

C++ + C — C + - — C_ + . . 

where the denominators in Eq.([3]) are the sum of all coincidences. In practice, the data 
{Ti, T 2 } are analyzed long after the data have been collected [8|. In general, the values for 
the coincidences C xy (a,f3) depend on the time-tag resolution r and the window W used to 
identify the coincidences, independent of which of the three pair identification procedures 
(see Section HVT) is being used. 



Data of EPRB experiments are often analyzed in terms of the function [a, 125J 



S{a, a', p, (3') = E(a, (5) - E{a, f3') + E{a', (3) + E(a', p'), (4) 

because it provides clear evidence that a quantum system is described by an entangled state. 
The idea behind this is that for any product state in quantum theory, or for the class of 
local realistic theories considered by Bell [uj 

- 2 < S(a, a', 13, (3') < 2, (5) 



an inequality known as one of Bell's generalized inequalities 25] . For later use, it is expedient 
to introduce the function 



S(6) = S(a, a + 29,a + 6,a + 39), 



(6) 



where we have fixed the relation between the angles (3 = a + 9, a' = a + 29, (3' = a + 39 
through the angle 9. Assuming rotational invariance, S(9) does not depend on a and we 
may set a = 0. 



III. QUANTUM THEORY 

In this section, we give a brief account of the quantum theoretical description of Cases I 
and II, strictly staying within the axiomatic framework that quantum theory provides. 

In the quantum theoretical description of Case I, the whole system is assumed to be 
described by the two-particle state 

= -j= (\HV) - \VH)) , (7) 

where H and V denote the horizontal and vertical polarization and the subscripts refer 
to photon 1 and 2, respectively The singlet state \^f) cannot be written as a product of 
single-photon states, hence it is an entangled state. 

In Case II, the photons have a definite polarization when they enter the observation 
station and the system is described by the product state 

|$) = (cos»7i|.H}i + sin 77i|F)i) (cos 772^)2 + sin r] 2 \V) 2 ). (8) 

The quantum theoretical expectation P+(a) (P_(/3)) for observing a photon at the + (— ) 
detector behind the polarizer with orientation a {(3) is given in the first two rows of Table [B 
The expressions for the two-particle correlation E(a,/3) are given in the third row. From 
Table [U it is clear that measuring E\{a) = P+(a) — P-(a), E 2 (P) = P+(P) — P-(fl) and 
E(a,(3) for various a and f3 suffices to distinguish between systems in the entangled state 
(Case I) or in the product state (Case II). 

In Case I, E(a, (3) — — cos 2(a — (3) and we find 

S{9) = 3 cos 29 -cos 69, (9) 
8 



TABLE I: The single- and two-particle expectation values for the two experiments described by 
the states Eqs. ([7]) and (jSJ), respectively. 

Case I Case II 

P+(a) 1/2 cos 2 (a-7?i) 

P-C9) 1/2 sin 2 (/?-r? 2 ) 

E\(a) cos2(a — 771) 

E 2 (f3) cos2(/3-r7 2 ) 

E(a,/3) — cos2(a — j3) cos 2(a — 771) cos 2(/3 — 772) 

which reaches its maximum value S max = maxg S(9) = 2\/2 at 9 = tc/8 + jir/2, where j is 
an integer number. 



Analysis of the experimental data , 110, 1 111 . Il2l |. yields results that are in 



good agreement with the expressions in Table HI leading to the conclusion that in a quantum 
theoretical description of Case I, the state does not factorize, in spite of the fact that the 
photons are spatially and temporally separated and do not interact. 

IV. DATA ANALYSIS OF A REAL EPRB EXPERIMENT WITH PHOTONS 

We analyze a data set (the archives Alice.zip and Bob. zip) of a real EPRB experiment 



with photons that is publicly available [26|]. The archives Alice.zip and Bob. zip contain 
data Eq. (pQ) for Case I for 9\ = 0,7r/4 and #2 = vr/8,37r/8. In a real experiment, such as 
the one described in Ref. the number of events detected at station 1 is unlikely to be 



the same as the number of events detected at station 2. The data sets of Ref. [26| show 
that station 1 (Alice.zip) recorded Aq = 388455 events while station 2 (Bob. zip) recorded 
A^ 2 = 302271 events. The fact that Aq 7^ A^ may have various reasons: It may happen that 
the source emitted one instead of two photons, one of the detectors did not respond to the 
arrival of a photon, the detector fired when there was no photon, etc. The data analysis 
does not account for such possibilities: We use the data as it is, without making additional 
hypotheses about unknown processes. 

We need a well-defined procedure to decide which two detection events stem from two 
particles that form a pair. Here, we use three different procedures to determine these 
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FIG. 3: (color online) Left: S max as a function of the time window W, computed from the data sets 



26l | . Squares 



contained in the archives Alice.zip and Bob. zip that can be downloaded from Ref. 
(black): Data obtained by comparing the detection times {i nj i|n = 1,...,N% = 388455} and 



{im^l'Ti = 1, . . . , N2 = 302271} with a reference clock. The maximum value of S„ 



2.78 is 



found at W = 4 ns at which the total number of coincidences is 2010 (~ 0.6%). Crosses (blue): 
Results of comparing the difference of the detection times with the time window W (A = 0). The 
maximum value of S max ~ 2.89 is found at W = 3 ns at which the total number of coincidences 
(with double counts removed) is 2899 (~ 0.8%). Bullets (red): Results of comparing the difference 
of the detection times with the time window W, taking into account the time shift A = 4 ns that 
maximizes the total number of coincidences, which (with double counts removed) is 13975 (~ 4%) 
in this case. The maximum value of S max ~ 2.73 is found at IF = 2 ns. Dashed line at 2\/2 ~ 2.82: 
Smax if the system is described by quantum theory (see Section J TTI ]) . Dashed line at 2: S max if the 
system is described by the class of models introduced by Bell [171 ] . Right: Same as left except for 
the range of W. 



pairs. Once the coincidences have been identified, we compute the two-particle average and 
Smax = maxgiS'(^) using Eqs.flH]), (JI]) and (EJ), respectively. In addition, we compute the 

frequency of coincidences, defined by 2(C ++ + C + C + _ + C_+)/ {N\ + jV 2 ). 

The first procedure divides the time interval of measurement (about 10 s) in equally 
spaced bins of size B For each station, we use the data {t n ^\n = 1, ...,JVi} and 

{tm,2\ m — 1) • • • > ^2} t° determine the number of events per bin and compute the coinci- 
dences by examining the content of the bins. This procedure compares the detection times 
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FIG. 4: (color online) Frequency of coincidences defined by 2(C ++ + C + C_| + C |_) / (iVi +-/V2), 

as a function of the time window VF (bin size B = 2W). The results were obtained by averaging 
the data of the four experiments (6>i = O,0 2 = tt/8), (6>i = 0,6*2 = 3tt/8), (0i = 7r/4, 6» 2 = vr/8), 



261 ] . Squares (black): 



and (6*i = ir/4,02 = 37r/8), contained in the archives Alice.zip and Bob. zip 
Data obtained by comparing the detection times {t n ^\n = l,...,Ni = 388455} and {t m ,2| m = 
l,...,iV2 = 302271} with a reference clock. Crosses (blue): Results of comparing the difference 
of the detection times with the time window W (A = 0). Bullets (red): Results of comparing 
the difference of the detection times with the time window W, taking into account the time shift 
A = 4 ns that maximizes the number of coincidences. For the values of W at which the three S max 
reach their maximum (squares: W = 4 ns, S max 2.78; crosses: W = 3 ns, S max ^ 2.89; bullets: 



W = 2 ns, S max 



2.73), the frequencies of coincidences are approximately 0.0015, 0.002, and 0.01, 



respectively. 



t n ^ and t m> 2 with a reference clock, using a coincidence window W = B/2. In the second 
procedure, we count the coincidences according to Eq. ([2]). In the third procedure, we ac- 
count for the fact that in the real EPRB experiment [9(], there may be an unknown shift A 
(assumed to be constant during the experiment) between the times £ n> i gathered at station 
1 and the times t m ^, recorded at station 2. Therefore, there is some extra ambiguity in 
matching the data of station 1 to the data of station 2. A simple data processing procedure 



11 



that resolves this ambiguity consists of two steps 9j. First, we make a histogram of the time 
differences t H) i —t m ^ with a small but reasonable time resolution (we used 0.5 ns). Then, we 
fix the value of the time-shift A by searching for the time difference for which the histogram 
reaches its maximum. Thus, we maximize the number of coincidences by a suitable choice 
of A. For the case at hand, we find A = 4 ns. 

The results for S max and the frequency of coincidences, as obtained by applying the three 
data analysis procedures, are presented in Figs. [3] and HI respectively. Note that in general, 
the frequency of coincidences depends on a and (3. However, for the choice a = 9\ = 0, 7r/4 
and (3 = # 2 = 7r/8, 37r/8, made in experiment by symmetry the four relevant frequencies 
of coincidences are expected to be the same, hence we show their average. As it is clear 
from Eq. (j2J) that the width of the time window in the second and third procedure is 2W, 
we have taken B = 2W to facilitate the comparison. 

From Fig. [3J, it follows that all three procedures yield a value of S mnx that significantly 
exceeds the upperbound (S max = 2) of the original Bell-like models [121]. As it has been 
shown rigorously that the original Bell (CHSH) inequality has to be modified if one uses 
Eq. (j2J) to count coincidences 141 ]. this violation should not come as a surprise. For W > 10 
ns and disregarding small fluctuations, the general trend is clear: S max decreases with W 
and drops below the "Bell-bound" for W > 300 ns. For W < 10 ns, each of the three 
procedures yields results for S max that are close to the quantum theoretical upperbound 
2^-2.83 |22|. 

The procedure that maximizes the coincidence count by varying A reduces the maximum 
value of S max from a value 2.89 (A = 0) that considerably exceeds the maximum (2y/2) for 



the quantum system [27|] to a value 2.73 (the value cited in Ref. [8j) that violates the Bell 



inequality 



171 ] and is less than the maximum for the quantum system. The fact that the 



"uncorrected" data (A = 0) violate the rigorous bound for the quantum system should not 
be taken as evidence that quantum theory is "wrong": It merely indicates that the way in 
which the data of the two stations has been grouped in two-particle events is not optimal. 

Analyzing the experimental data set 26j with A = 4ns and W = 2ns (yielding the "best" 
value of S max ~ 2.73 and a total number of coincidences of 13975) gives £i(0, n/8) = —0.073, 
£i(0,3vr/8) = 0.118, £i(vr/4, tt/8) = 0.036, £i(tt/4, 3tt/8) = -0.065, and E 2 (0,vr/8) = 
0.188, £ 2 (0,37r/8) = 0.258, £ 2 (tt/4, tt/8) = 0.099, £ 2 (tt/4, 3tt/8) = -0.147, significantly dif- 
ferent from the theoretically expected value (zero, see Tabled]). Disregarding the coincidence 
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criterion, we find J2 n x n,i = -0.007,-0.005 for a = 0, vr/4 and ^ n i n , 2 = -0.028,-0.024 
for (3 = 7r/8,37r/4. Apparently, all these numbers change considerably with the station and 
with the settings of the electro-optic modulators. 

The results presented in Fig. HI show that for small W, the frequency of coincidences 
depends significantly on the criterion that is used to identify pairs of events. The procedure 
that maximizes the coincidence count (procedure three) seems to yield the most "stable" 
results. For all three procedures, the frequency of coincidences at which S max reaches its 
maximum is below 1% (see Fig. HJ). If we identify the observed frequency of coincidences 
with the probability of coincidences 7_that enters the upperbound in the properly modified 
Bell inequality (see Eq. (16) in Ref. [14|), then this theoretical upperbound is larger than 4, 
supporting the idea that the experimental data {§1 is not in conflict with local realism [^J. 

Textbook treatments of EPRB experiments assume that the correlation, as measured in 
experiment, is given by [13] 

N 

c if = ( 10 ) 

n=l 

which is obtained from Eq. (j2J) by taking the limit W — > oo, hence the notation Ciy. 
Although the limit W — > oo defines a valid theoretical model, there is no reason why this 
model should have any bearing on the real EPRB experiments with photons as they have 
been performed so far. An argument that might justify taking the limit W — > oo is the 
hypothesis that for ideal experiments, the value of W should not matter. However, as our 
analysis of the experimental data shows, to make contact to quantum theory, one has to 
reduce (not increase) W js, 9]. Thus, in real EPRB experiments with photons, the window W 



matters (a,|9j. The details of the criterion that is used to decide which two events correspond 
to the observation of a single two-particle system seem to be of secondary importance. 

As it is relatively easy to reproduce the results of quantum theory in the regime of small 
W [ijj], and as keeping W arbitrary does not render the mathematics more complicated, 
there really is no point of studying the simplified model defined by Eq. (jlOp : We may always 
consider the limiting case W — > oo afterwards. 
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V. SIMULATION MODEL 



We now take up the main challenge, the construction of Einstein-local, causal processes 
that generate the data sets Eq. (CQ) such that they reproduce the results of quantum theory, 
summarized in Table [H A concrete simulation model of the EPRB experiments sketched 
in Figs. [H and [2] requires a specification of the information carried by the particles, of the 
algorithm that simulates the source and the observation stations, and of the procedure to 
analyze the data. From the specification of the algorithm, it will be clear that it complies 
with Einstein's criterion of local causality on the ontological level: Once the particles leave 
the source, an action at observation station 1 (2) can, in no way, have a causal effect on the 
outcome of the measurement at observation station 2 (1). 



A. Source and particles 

The source emits particles that carry a vector S n)i = (cos(£ n + (i — l)7r/2), sin(£ n + {i — 
l)7r/2)), representing the polarization of the photons. The "polarization state" of a particle 
is completely characterized by £ n , which is distributed uniformly over the whole interval 
[0, 2n[. We use uniform pseudo-random numbers to mimic the apparent unpredictability of 
the experimental data. However, from the description of the algorithm, it trivially follows 
that instead of uniform pseudo-random number generators, simple counters that sample the 
intervals [0, 2ir[ in a systematic, but uniform, manner might be employed as well. This is 
akin to performing integrals by the trapezium rule instead of by Monte Carlo sampling. The 
source thus emits two particles with mutually orthogonal, random polarization. 

In Case II we change the unpredictable polarization state of the particles to a fixed polar- 
ization state by placing polarizers in between the source and each observation station. These 
polarizers have one input and one output channel and their orientations are characterized 
by the angles rji and rj2. 



B. Observation stations 



Prior to collecting data, we fix the number M of different polarization directions (M = 2 
in the experiment of Ref. [8|). We use 2M random numbers to fill the arrays («i, «m) and 
(3m)- Before (or after) the nth pair leaves the source, we use two uniform random 
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numbers 1 < m,m' < M to select the angles 6 n i = a m and 9 n ^ = Pm'- In practice, we use 
two different pseudo-random number generators for observation stations 1 and 2, but we 
have never seen any statistically significant effect of using the same one for both observation 
stations. 



C. Polarizer 



We make the hypothesis that in laboratory EPRB experiments with photons the various 
polarizers are interchangeable. Therefore, the algorithm to simulate the two polarizers in 
the observation stations should be identical. Evidently, for the present purpose, if we switch 
from Case I to Case II, it is not permitted to change the algorithm for the polarizer. This 
also holds for the polarizers placed in between the source and the observation stations. 

The input-output relation of a polarizer is rather simple: For each input event, the 
algorithm maps the input vector S onto a single output bit x. The value of the output bit 
depends on the orientation of the polarizer a = (cos a, sin a). According to Malus law, for 
fixed S = (cos £, sin £) and fixed a, the bits x n are to be generated such that 

1 N 

i im T?y2 Xn = cos2 (£ - «)> ( n ) 

n=l 

with probability one. If the input vectors S are distributed uniformly over the unit circle, 
the sequence of output bits should satisfy 

1 N 

71=1 

with probability one, independent of the orientation a of the polarizer. 
The model for a polarizer is defined by the rule 

{+1 if r n < cos 2 (£„ - a) 
-1 if r n > cos 2 (£ n - a) 

where < r n < 1 are uniform pseudo-random numbers. The polarizer sends a particle with 
polarization S n j = (cos a, sin a) or S nj j = (— sin a, cos a) through its output channel +1 or 
— 1, respectively. It is easy to see that for fixed £ n = £ and a, this algorithm generates events 
such that 2(x n> i) — 1 = cos 2 (£ — a), where (X) denotes the average over many realizations 
of the variables r n and £ n . In this case, the input-output relation of the simulation model 
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agrees with Malus law Eq. (II ip . On the other hand, if £ n is distributed uniformly over the 
interval [0, 2ir[, we have (x n> i) = 0, in agreement with Eq. (fT2"j) . It is at this point, the model 



for the polarizer, that the simulation model differs from the one used in Ref. H5|: The model 
of the polarizer used in Ref. [15| can reproduce the correlation of the singlet state but cannot 
reproduce Malus law. 

In Case II we discard particles with polarization r/i + tt/2 (r] 2 + vt/2) that leave the 
polarizers placed in between the source and observation station 1 (2). 



D. Time delay 

In our model, the time delay t n ^ for a particle is assumed to be distributed uniformly over 
the interval [to,t + T]. In practice, we use uniform pseudo-random numbers to generate 
t n ^. As in the case of the angles £ n , the random choice of t n ^ is merely convenient, not 
essential. From Eq.([2]), it follows that only differences of time delays matter. Hence, we 
may put to = 0. The time-tag for the event n is then t n ^ e [0, T]. 

There are not many reasonable options to choose the functional dependence of T. As- 
suming that the particle "knows" its own direction and that of the polarizer only, T should 
be a function of the relative angle only. Furthermore, consistency with classical electrody- 
namics requires that functions that depend on the polarization have period ir [28J. Thus, 
we must have T(£ n — #1) = F((S nj i • a) 2 ) and, similarly, T(£„ — 9 2 ) = F((S nj 2 ■ b) 2 ), where 
b = (cos /3, sin (3). We found that T(x) = T |sin2x| d yields the desired results [15(. Here, 
T = maxfl T(6) is the maximum time delay and defines the unit of time, used in the simu- 
lation. In our numerical work, we set T = 1. 



E. Data analysis 

For fixed N, the algorithm described earlier generates the data sets T i; just as experiment 
does. In order to count the coincidences, we choose a time-tag resolution < r < T and 
a coincidence window W > r. We clear all the coincidence counts C xy (a m , (3 m i) for all 
x, y — ±1 and m,m' = 1,...,M. Then, we make a loop over all events. To count the 
coincidences, we first compute the discretized time tags k nti = [t„,i/r] for all events in 
both data sets. Here \x] denotes the smallest integer that is larger or equal to x, that is 
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\x\ — 1 < x < \x\. According to the procedure adopted in the experiment [8(, an entangled 
photon pair is observed if and only if |fc n l — k n>2 \ < k = \W/t}. Thus, if \k njl — k n>2 \ < k, 
we increment the count C XnAiXn2 (a m , (3 m >). 

We emphasize that the simulation procedure counts all events that, according to the same 
criterion as the one employed in experiment, correspond to the detection of two-particle 
systems. Note that in our simulation model, the three different methods that we used to 
analyze the experimental data (see Section IPVT) give identical results. 



VI. PROBABILISTIC TREATMENT 

Let us assume that we can analyze our simulation model, described in Section |V] by 
replacing the deterministic sequence of pseudo-random numbers by the mathematical con- 



cept o 



ity 
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independent random variables, as defined in the (Kolmogorov) theory of probabil- 



30] . Under this assumption, each event constitutes a Bernouilli trial [29j, 130J and we 
can readily obtain analytical expressions for the expectation values that we compute with 
the simulation model. 

This section serves three purposes. First, it provides a rigorous proof that for up to first 
order in W and for d — 4, the probabilistic description of the simulation model exactly 
reproduces the single particle averages and the two-particle correlations of quantum theory 
for the system under consideration. Second, it illustrates how the presence of the time- 
window introduces correlations that cannot be described by the original Bell-like "hidden- 
variable" models [yj]. Third, it reveals a few hidden assumptions that are implicit in the 
derivation of the specific, factorized form of the two-particle correlation that is essential to 
Bell's work. 

As explained in Section [Til rea l EPRB experiments with photons produce the data sets 

Ti = {x nt i = ±1, t nti , 9n ti \n = 1, . . . , N} . (14) 

Let us assume that there exists a probability, denoted by P(xi,x 2 ,ti,t 2 \a, j3), to observe 
the data {xi,ti, 0\ = a} and {x 2 , t 2 ,9 2 = ft} at station 1 and 2, respectively. Notice that we 
assume, unlike in the computer simulation model where 9 n ^ may change with each event n 
but as in the case of quantum theory, that a and /3 are fixed. The mathematical expectation 
of the coincidences C xy (see Eq. ([2])), that is the average computed with P(xx, x 2 , ti, t 2 \a, 
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is given by 



/+oo c+oo 
dh dt 2 P(x,y,t 1 ,t 2 \a,P)Q(W- |*i-t 2 |). (15) 
-oo J — oo 

Once we know (C xy ), the mathematical expectation of the single-particle counts and two- 
particle coincidences follow from 



E 2 (a,(3,W) 



Ei,y=±l X (C X y) 

y{C X y) 

Y2x,y=±l(C X y) 
Y^x,y=±l X V{Cxy) 



E(a,(3,W) = ^=- r (16) 

As a first step, let us express the probability for observing the data {xi, x 2 , tj, t 2 } as 
an inte gra l over the mutually exclusive events £i,£ 2 - According to the rules of probability 



theory [29|, |30|] , we have 



2tt r 2n 



P(x 1 ,x 2 ,t u t 2 \a,{3) = J J P{x l ,x 2 ,t l ,1*\a i p i Z l ,&)P(Z l ,&\a,P)d£ l d£2 1 (17) 

where £i and £ 2 denote the two-dimensional unit vectors, representing the polarization. 
Starting from the exact representation Eq. (1171) . we now assume that in the probabilistic 
version of our simulation model, for each event, the values of {xi,x 2 , ti,t 2 } are independent 
of each other and that the values of {xi,ti} ({x 2 ,t 2 }) are also independent of f3 and £2 (a 
and £1)). Thus, we may write 

I /*2"7r /•27T 

P(x 1 ,x 2 ,t 1 ,t 2 \a,p) = — J J P(x 1 ,t 1 \x 2 ,t 2 ,a,f3,£ 1 ,£ 2 )P(x 2 ,t 2 \a,p,£ 1 ,£ 2 ) 

xP(£i,£ 2 \a,(3)d£id£ 2 
P{x u h\a, £i)P(x 2 , t 2 \(3, e 2 )P(6, £ 2 |«, (3)d^ 2 



1 r 2w r2ir 



4vT 2 
1 

1 

4^ 



Jo 

2vr p2ir 



JO 



P( Xl \a, Zi)P{h\a, £i)P{x 2 \P, &)P(t 2 \P, £ 2 ) 

*P(£i,&\a,P)d£id£ 2 



P{xi\a, £i)P(*i|a, £i)P(x 2 |/3, £ 2 )P(* 2 |/3, £2) 



Jo 



where, in the last step, we assumed that the values of £1 and £ 2 are independent of a or /3. 
With the three assumptions made so far, Eq. (|18|) gives the exact probabilistic description of 
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our simulation model. It is of interest to note that Eq. (lisp can be derived directly from the 
description of the algorithm, without recourse to probability theory, by letting the number 
of events in the discrete sums approach infinity. 

The mathematical structure of Eq. ( fl8l) is the same as the one that is used in the derivation 
of Bell's results and if we would go ahead in the same way, our model also cannot produce 
the correlation of the singlet state. However, the real factual situation in the experiment 8| 
is different: The events are selected using a time window W that the experimenters try to 
make as small as possible Q| . Accounting for the time window, that is multiplying Eq. (flHj) 
by the step function and integrating over all t x and t 2 , the expression for the probability for 
observing the event (x x ,x 2 ) reads 

f C P{xi\<*, £i)P(x 2 \P, £ 2 )w(a, P, 6, 6, W)P(£ U 6)^1^2 



P(x x ,x 2 \a,f3) 



(19) 



where, in general, the weight function 

/+00 r+00 
dt x / dt 2 P{t x \a, £i)P(h\P, 6)Q(^-I*i-*2|), (20) 
-00 ./ —00 

will be less than one (because fj 1 ^ dti J^°° dt 2 P(t x \a, £ x )P(t 2 \f3, £2) = 1) unless W is larger 
than the range of (t x ,t 2 ) for which P(t x \a,£ x ) and P(t 2 \f3,!; 2 ) are nonzero. It is self-evident 
that unless w(a, (3, £1, £ 2 , W) = w(a, £ x ,W)w((3,£ 2 ,W), Eq. (JTHJ) cannot 
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De written in the 
for the nota- 
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factorized form P(xi, x 2 \a, j3) = f P(x l \a, X)P(x 2 \f3, X)p(X)d\ (see Ref. 
tion) that is essential to derive the original Bell (CHSH) inequalities. 

In our simulation model, the time delays are distributed uniformly over the interval 
[0,Ti] where T\ = T |sin2(a - £i)| d and T 2 = T |sin2(/? - £ 2 )| d - Thus, P(t x \a,£ x ) = 

e(t x )e(T x -t x )/T h P(h\P,&) = e(t 2 )e(r 2 - t 2 )/r 2 , and 

w(a,p,^ 2i W) = r dt x f 2 dt 2 ®{W-\t x -t 2 \). (21) 

Jo Jo 

The integrals in Eq. fl2"TT) can be worked out analytically, yielding 

w(a, (3, Ci, 6, W) = [ Tl + 7| + 2(T X + T 2 )W + (W - T x ) \W - T x \ 

*i± X ± 2 

+(W -T 2 )\W-T 2 \-(W-T X + T 2 )\W-T X +T 2 \ 
-(W + T X -T 2 )\W + T X -T 2 \ ]. (22) 
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Clearly, Eq. (I2"2"j) cannot be written in the factorized form w(a, £1, W)w(fl, £ 2 , W 7 ). Hence, it 
should not come as a surprise that as soon as we want to simulate the real EPRB experiment 
with photons in which the time window is essential, we can obtain correlations that cannot 
be described by Bell-like models. 

According to our simulation model (and the assumption made at the beginning of this 
section), the probability distributions that describe the polarizers are given by 

1 + xi cos 2(a — £1) 



2 



tdi \at\ 1 + ^2 cos2(/?-£ 2 ) f0Q s 
P(x 2 \f3,&) = • (23) 

It is easy to check that these distributions reproduce Malus law for a single polarizer. 

We now consider some specific cases. First, we consider Case I and specialize to the case 
that the source emits particles with opposite polarization P(£i,£ 2 ) — + 7r /2 — ^)-P(^i) 
with P(£i) being a uniform distribution. If d = and W < T , we have w(a, /3, £i, £ 2 , W) = 
(2T - W)W/Tl Likewise, if W > T , w(a, (3, £ 2 , W) = I. Therefore, if W > T or d = 0, 
we have 

P{X\,X<2\Oi.p) = — 7, 7, 

P( Xl \a, ii)P{x 2 \(3, &)P{1- U 6)^1^2 



JO 



2k 



1 + x x cos 2(a - - «2 cos 2 (/3 - £))*; 



1 

8^ 

2 + x\X2 cos 2 (a — (3) 



8 



(24) 



showing that if we ignore the time-tag information, the two-particle probability takes the 
form of the hidden variable models considered by Bell 17|, and we cannot reproduce the 



results of quantum theory 



Second, we consider Case I but focus on the regime for small W, the regime that experi- 
menters aim to reach j^. Then, Eq. ( !22l) reduces to 



2W 

w(a, 0, fc, 6) = TTfr^T + OiW 2 ). (25) 

max(Ti, T 2 ) 
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and we find that Ei(a,(3) = E 2 (a,(3) = and that 

J 2 ^cos2(£ - a) cos2(£ - 0) max(| sin2(£ - a)|, | sin2(£ - P)\y d d£ 



E(a,/3) 



f maxfl sin 2 - a) |, | sin2(£ - f3)\)- d d£ 
/; / / 2 2+7r/4 cos2ecos2(e - 9)\ sin2^g 

;;/ 2 2+ ^lsin2ei-^ 

' -±cos2fl ,d = 

I sin 29 cos 29 -cos 29 + In [| tan 9 \ si ^ 2e / 2 } } d=2 

- cos 29 , d — 4 , 
-i cos 20 [1 + 24(19 + 5 COS40)- 1 ] ,d=6 

- (53 cos 29 + 7 cos 69) (39 + 21 cos A9)~ l ,d = 8 



(26) 



where 9 = a — (3 and we have omitted the expressions for odd d because they cannot be 
written in terms of elementary functions. Note that by passing to the continuum limit, one 
has to be careful with the integrals that appear in expressions such as Eq. (1261) : The ratio of 
the two integrals is well-defined but each individual integral may vanish or diverge for some 
choices of a — (3. Needless to say, such situations do not occur in the (discrete) simulation 
model. 

Finally, in Case II, the values of £i and £2 are fixed, hence P(£i, £2) = <K£i — Vi)3(£,2 ~ ^2)- 
Then, as is clear from Eq. (fl9l) . the weight function w(a, /3, rji, 772, W) drops out and the 
two-particle probability reduces to 



P(xi,x 2 \a,P) = P(xi\a,T)i)P(x2\P,7}2), 



(27) 



such that 



E(a,{3,W) = ^2 x 1 x 2 P(xi,x 2 \a, P) = ^ x 1 P(x 1 \a, (3) \ ^ x 2 P(x 2 \a, (3) 

Xl,X 2 = ±l \Xl=±l / \X 2 = ±1 / 



E l (a,(3,W)E 2 {a,f3,W). 



(28) 



Evidently, the simulation model will reproduce the results of quantum theory for Case II if 
the proper expression, the one yielding Malus law, is used for the single-particle probabilities 
P(x 1 \a 7 rji) and P(x 2 \/3 J r] 2 ). 

Summarizing: Up to first order in the time window W and for d = 4, in Case I (corre- 
sponding to the case in which the source emits particles with opposite random polarization) 
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the probabilistic model of the simulation algorithm yields 



E 1 (a,/3) = E 2 (a,P) = 

E{a,/3) = -cos2(a-/3), (29) 

for the single-particle averages and two-particle correlation, respectively. Obviously, these 
expressions are identical to those given in the second column of Table [B If, as in Case II, 
the source emits particles with fixed polarizations 771 and r/2, respectively, the probabilistic 
model of the simulation algorithm yields 

Ei(a,P) = cos2(a — 771) , E 2 (a, (3) = cos2(a — r/ 2 ) 
E(a,/3) = — cos2(a — 77x) cos2(a — ?7 2 ), (30) 

in exact agreement with the results in the third column of Table [H Thus, it follows that to 
first order in W, the probabilistic model of the simulation algorithm can reproduce exactly 
the results for the single- and two-particle averages of the quantum theory of a system of 
two photon polarizations. 

VII. SIMULATION RESULTS 

We use the computer model, described earlier to simulate Cases I and II. The simulation 
proceeds in exactly the same way as in the experiment, that is we first collect the data 
sets Ti and T 2 for various settings of the polarizers (various 9 nji ), and then compute the 
coincidences Eq.((2]), the average single-particle counts and the correlation Eq.Q, from which 
we can calculate the function S(6) (see Eqs. (HI) and ([61)). The parameters for all simulations 
are k = 1, d — 4, r = 0.00025, and N = 10 6 , unless mentioned otherwise. 

In Fig. [5] (left), we present simulation data for the correlation E(ct, (5) for Case I, that is for 
the case that the source emits particles with an opposite, random polarization, corresponding 
to the singlet state in the quantum theoretical description. Figure E] (right), shows the 
corresponding results for Case II. It is clear that in both cases, the agreement between the 
simulation data and quantum theory is excellent. 

Also shown in Fig.Oare the results for E(a, (3) if we ignore the time-delay data (equivalent 
to d = or W — * 00). In Case I we obtain simulation results that agree very well with the 
expression E(a,/3) = —(1/2) cos 28 (see Eq.f l26p ). a result that differs from what is obtained 
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FIG. 5: (color online) Correlation E(a, f3) between the coincidence counts as a function of the 
orientation difference of the two polarizers in each observation station. Left: Computer simulation 
of Case I in which the source emits particles with opposite random polarization (EPRB experiment). 
Right: Computer simulation of Case II in which the source emits particles with fixed polarization 
and a = 9, = 9 + 7r/4, 771 = ir/6, and 772 = 7r/6 + ir/2 (see Fig. [2]). Squares (red): Simulation 
results using the time-delay mechanism (with d = 4) to compute the two-particle coincidence. Open 
circles (black): Simulation results without using the time-tags (equivalent to d = or W — > 00). 
Solid lines: Quantum theory. 



by considering the class of models studied by Bell [13]. In the latter case an equilateral 
saw-tooth function is obtained instead of the cosine. In Case II, the results for d — 4 and 
d = or W — > 00 are, apart from statistical fluctuations, the same. Hence, for Case II the 
time window W can be omitted for the calculation of the two-particle correlation function. 

In Fig. [U] (left), we present additional simulation data for Case I. It is clear that for 
d = 4, the simulation model reproduces the results of quantum theory for the single-particle 
expectation values P±(a) and P±(/3) (see Table fl} and S(9) (see Eq. ([9])). Indeed, the 
frequency with which each detector fires is approximately one-half and |5'(^)| agrees with 
the expressions E(a,/3) = — cos2(a — 0) that is obtained for the singlet state. Also shown 
in Fig. E] (left) are the results for |if>(0)| if we ignore the time-tag data. Effectively, this is 
the same as letting the time window W — > 00 or setting d = 0. Then, our simulation model 
generates data that satisfies | S(9)\ < 2, which is what we expect for the class of models 
studied by Bell [17]. 

In Case II, the source emits particles with a fixed (but not necessarily opposite) po- 
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FIG. 6: (color online) Left: Computer simulation of Case I in which the source emits particles 
with opposite random polarization (EPRB experiment). Right: Computer simulation of Case II 
in which the source emits particles with fixed polarization. Squares (red): Simulation results for 
F{9) = \S{9)\ using the time-delay mechanism (d = 4). Open circles (black): Simulation results 
for F{9) = | S{9)\ without using the time-tags (equivalent to d = or W — > oo). Other markers: 
Average single-particle counts on the detectors (see Fig. Q]). Squares (green): F(9 = 9i) = P + (9i); 
Diamonds (green): F{9 = 9 X ) = P_(0i); Plusses (blue): F{6 = 9 2 ) = P+(9 2 ); Crosses (blue): 
F{9 = #2) = P-{92). In Case I (left), these four symbols lie on top of each other. In Case 
II (right), these markers show the typical Malus law behavior. Solid line: Quantum theory for 
| S(9) | . Dashed line at \S(9)\ = 2y/2: Maximum of S(9) if the system is described by quantum 
theory. Dashed line at | <!?(#) | = 2: Maximum of S{9) if the system is described by the class of 



models introduced by Bell 17]]; Dashed line at |5'(0)| = 1/2: Expected number of +1 and — 1 
events recorded by the detectors if the input to the polarizers consists of particles with random 
polarization. Dotted lines: Quantum theory for P + {9\), P_(9i), P + {92) and P_(#2)- 



larization. In the right panel of Fig. El we present results for the case 9\ = a = 9 and 
#2 — ft = 6-\-7i/4. The angle £ of the particles is 7r/6 (corresponding to rji = 71/6 
and 772 = 7r/6 + 7r/2 in the quantum theoretical description). For this choice, we have 
P + (a) = cos 2 (0 - tt/6), P+(ft) = cos 2 (6 - tt/6 - tt/4), E{a,ft) = 2- 1 sin4(7r/6 - 9) and 
S{6) = sin4(7r/6 — 6). Also seen from Fig. El (right) is that |5'(6')| does not depend on d, 
that is apart from statistical fluctuations, the time-tag data do not affect 

From Fig. El it is clear that for d = 4, the event-by-event simulation model reproduces 
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FIG. 7: (color online) Left: Same as Fig. © (left) except that d = 2. Dotted line: F(0) = \S(9)\ 
calculated from Eq.{2B]). Right: Same as Fig. © (left) except that d = 6. Dotted line: F(0) = \S(6)\ 
calculated from Eq. (|26l) . 

the single- and two-particle results of quantum theory for both Case I and II, without any 
change to the algorithm that simulates the polarizers. 

Having established that the data generated by our "non-quantum" system agrees with 
quantum theory, it is of interest to explore if these algorithms can generate data that is 
not described by quantum theory and by the locally causal, probabilistic models introduced 
by Bell [jjj. We can readily give an affirmative answer to this question by repeating the 
simulations for Case I (see Fig. [6] (left)) for different values of the time-delay parameter d, 
all other parameters being the same as those used to obtain the data presented in Fig. [6] 

For d = 0, simulations with or without time-delay mechanism yield data that, within 
the usual statistical fluctuations, are the same (results not shown) and satisfy | S(9)\ < 2 . 
Figure [7] shows the simulation data for d = 2 and d = 6. For < d < 4 our model yields two- 
particle correlations that are stronger than those of the Bell-type models but they are weaker 
than in the case of the singlet state in quantum theory. Therefore, the maximum of S(6) is 
less than 2y/2 but larger than two. For d > 5, we find that the two-particle correlations are 
significantly stronger than in the case of the singlet state in quantum theory. 

For d < 4, 2 < S max < 2y/2 for any value of W/r. Hence, for d < 4 our model 
cannot produce the correlations of the singlet state. For d — 4, 2 < S max < 2\/2 and 
our model produces the correlations of the singlet state if W is sufficiently small such that 
contributions of order W 2 can be negelected. For d > 4, 2 < S max < 4, and for a range of 
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FIG. 8: (color online) Left: S max = max^Q,/^/?/ S(a, a', (3, /?') as a function of the time window 
W relative to the maximum time delay resolution To. Curves from bottom to top: Results for 
d = 0,1,..., 10. Solid line at 2y/2: Value of S max for a quantum system in the singlet state. 
Dashed line at 2: Value of S max for a quantum system in an uncorrelated state. Right: S max = 
max a a i g gi S(a,a' , /3, (3') as a function of the time-delay parameter d negelecting contributions of 
OiW). Solid line at 2\/2: Value of S max for a quantum system in the singlet state. Dashed line at 
2: Value of S max for a quantum system in an uncorrelated state. 

W/t, S max > 2y / 2, implying that our model exhibits correlations that cannot be described 
by the quantum theory of two spin-1/2 particles, while still satisfying Einstein's criteria for 
local causality. 

From Fig. [7J it can be seen that for d = 2 and d = 6 there is good agreement between the 
results obtained with our event-based simulation model and the analytical result for |if>(0)| 
obtained from Eq. (1261) . For d = 2, the simulation results show larger fluctuations than for 
d = 6, but in all cases these fluctuations can be reduced by increasing N (results not shown). 

The simulation results presented in Figs. E] and [7] have been obtained for W/t = 1 and 
small t (recall that the unit of time in our numerical work is set equal to one). In general, 
in experiment the two-particle correlation depends on both W and r. Our simulation model 
makes definite predictions for this dependence. This can be seen from Fig. [HI showing 
Smax = maXa,a',/),/3' S(a, a', (3, (3') as a function of W/T for various values of d and of S max 
to first order in W as a function of d. The numerical results agree with the values of S max 
that have been obtained analytically to first order in W, d = 0, 2, 8 and for W — > oo. 

Summarizing: In the regime of small W/T = r/T , the results produced by the simulation 
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algorithm are in excellent agreement with the quantum theoretical expressions (see Table [T|) 
of the single- and two-particle averages. 



VIII. COMPARISON WITH EXPERIMENTAL DATA 

For the simulation model described in Section IVJ it follows immediately from Eqs. (J2]) 
and (131) that \E(a,j3)\ < 1 and that 

\E(a, (3) - E(a, (3') + E(a', (3) + E(a', (3') | < 4. (31) 

Without any further constraints on the algorithm that generates the data {Y 1 ,Y 2 } the 
upperbound (4) in Eq. f[3"Tj) cannot be improved. On the other hand, for a local realist 
(probabilistic) model, it can be shown that [ljj] 

\E(a, (3) - E(a, (3') + E(a', (3) + E(a', /3') | < - - 4, (32) 

7 

where 7 is the infimum of the probability of coincidence over all possible settings {a, j3}. In 
our simulation model the frequency of coincidences 

1 - 

r = -^ e> ( M/ -IW-Wl) ) (33) 

n=l 

is easy to compute and, assuming that the results that we obtain by using pseudo-random 
numbers can be described by a probabilistic model (see Section IVT1) . we may assume that 
7 = T with probability one. For W = t and d = 4, a straightforward calculation gives 

7 = minr = ^ + 0(^ 2 ), (34) 
q-/3 on 1q 

showing that up to first order in the time window W , the minimum frequency of coincidences 
is proportial to the time window, as one naively would expect. 

As pointed out in Ref. [iuj . a local realist model that uses coincidence in time to decide 
which particles form a pair is not necessarily in conflict with the predictions of quantum 
theory unless 7 > 3 — 3/v2. Thus, in Case I, it is of interest to explore how T affects S max 
and the sinusoidal shape of the two-particle correlation but before we present some results, 
we want to draw attention to the fact that the model that we introduce in this paper is not 
unique in the sense that it is not the only model that reproduces the results of quantum 
theory for the singlet state [ijj]. Different models will yield different numerical results for 
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r but the general behavior is the same. In these event-based models, there are three inde- 
pendent parameters that we can use to "tune" the simulation results to experimental data, 
namely r/T , W/T , and d. 

First, we consider the problem of "fitting" our model to experimental data for S max . 
As explained earlier, to obtain S max , one has to perform four experiments. For instance, 
Ref. [8] reports S max ~ 2.73 using a = Q\ = 0,7r/4 and ft = 9 2 = 7r/8,37r/8 and if we 
make the hypothesis that the frequency of coincidences that we found earlier 0.01 for 
a — ft = ±7r/8, — 37r/8) is a good estimate for 7, no conclusion can be drawn from the 
relevant theoretical bound Eq. ( |32|) . other than that this experiment does not rule out a 
local realist (probabilistic) description [jjj] . The simulation model described in Section IVl 



reproduces the experimental result S max ~ 2.73 [8( for r/T = 1/29, W/T = 1/29, and 
d = 4. For a -ft = 0, tt/8, tt/4, 3tt/8, tt/2 we find T = 0.38, 0.14, 0.06, 0.14, 0.38, respectively. 
Because in our simulation the source only emits particle pairs and since no particles are lost 
or falsely detected, we may expect that the simulation for a — ft — ±7r/8, — 37r/8 yields a 
value of T that is larger than the one (~ 0.01) extracted from the experimental data jsj]. 
A simulation run with N=300000 events (roughly the same number as observed in the 
experimental data analyzed in Section HVT) . gives J2 n x n,i = 0.0016, —0.0011 for a = 0,7r/4 
and Ylm x n,2 = 0.007,0.001 for ft = 7r/8,37r/4, respectively, in reasonable agreement with 
the experimental results (see Section [IVl) . 



For comparison, the simulation model introduced in Ref. 15|] reproduces the same value 
of S max for r/T = 1/9, W/T = 1/9, and d = 2, yielding T = 0.38,0.18,0.13,0.18,0.38 for 
a — ft = 0, 7r/8, 7r/4, 37r/8, tt/2, respectively. 

As another example, we consider the result S max ~ 2.25 as obtained from ion-trap ex- 
periments 1Q( |. Although it is not evident that the events registered in this experiment 
are as simple as the detection of single photons, let us assume that the model for the 
real EPRB experiment with photons can nevertheless be used to describe the outcome 
of these ion-trap experiments. Then, the simulation model described in Section |V] repro- 
duces the value of S max ~ 2.25 if we take r/T = 1/4.3, W/T = 1/4.3, and d = 4. For 
a - ft = 0, tt/8, tt/4, 3tt/8, tt/2 we find T = 0.65, 0.46, 0.32, 0.46, 0.65, respectively. For com- 
parison, the simulation model described in Ref. [15| yields the same value of S max = 2.25 
for r/T = 1/1.031, W/T = 1/1.031, and d = 2, with T = 0.95,0.89,0.89,0.89,0.95 for 
a — ft = 0, 7r/8, 7r/4, 37r/8, 7r/2, respectively. As T w 1, this experiment seems to have an 
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FIG. 9: (color online) Left: The frequency of coincidences r as a function of \a — (3\ for parameters 
t/Tq, W/Tq and d chosen such (see text) that the simulation model reproduces quantum theory 
(solid line, red), S max = 2.83, and the values of S max = 2.25 (dashed line, black) and S max = 2.73 
(dotted line, blue), as obtained from experiments with ions K] and with photons [8J, respectively. 
Right: Simulation results of the two-particle correlation E(a — (5) as a function of \a — /9| for 
the model parameters that yield S max = 2.25 (squares, black) and S max = 2.73 (stars, blue), 
respectively. The dashed lines are given by —0.875 cos 2(a — /?) and —0.725 cos 2(a — (3). The solid 
line (red) is the result — cos2(a — j3) of quantum theory (see Fig. [5] (left)). 



almost ideal detection efficiency [10 ]. 

For completeness, we consider the case S max = 2.83. Recall that both the model intro- 
duced in this paper and the one of Ref. [15J reproduce the result (E(a — (3) = — cos2(a — (3)) 
of quantum theory if we keep the contributions of 0(W) only. The model described in Sec- 
tion M yields S max = 2.83 if we take r/T = 1/1500, W/T = 1/1500, and d = 4 for which 
r = 0.13,0.0032,0.0011,0.0032,0.13 for a - (5 = 0, tt/8, tt/4, 3tt/8, tt/2, respectively. For 
these choices of parameters, the numerical results for E(a — (3) are very close to those of 
quantum theory (see Fig. El). 

For the model of Ref. [15J and r/T = 1/1500, W/T = 1/1500, and d = 2, we find 
r = 0.031, 0.0011, 0.00085, 0.0011, 0.031 for a - [3 = 0, tt/8, tt/4, 3tt/8, tt/2, respectively. 

In Fig. [9]4eft), we plot T as a function of a — (3, as obtained for the simulation model 
introduced in the paper, for the three cases S max = 2.25, 2.73, 2.83 discussed earlier. The 
general trend is clear: T reaches its maximum at a — (3 = 0, 7r/2, ... and its (nonzero) 
minimum at a — (3 = 7r/4, 37r/2, . . .. 
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Finally, we study how E(a — (3) deviates from the result E(a — /?) = — cos2(a — (3) of 
a system in the singlet state as we fit the values of S max to the experimental results. In 
Fig. fright), we show the simulation results for the two cases S max = 2.25, 2.73, correspond- 
ing to the experiment with ions of Ref. Q and the experiment with photons of Ref |8|], 
respectively. From Fig. fright), we see that the main effect of reducing S max is to reduce 
the amplitude (visibility) of the correlation. Although it is clear that the simulation data 
cannot be described by a single sinusoidal function, the deviations are small and it remains 
to be seen if experiments can resolve such small differences. 

As is evident from Fig. [BJ for d > 4 our model yields the value for the singlet state 
Smax = 2 a/2 without having to consider the regime of small W. Thus, in order for an 
experiment and a model of the type considered in our paper to reproduce the features of 
a quantum system of two 5=1/2 particles in the singlet state, it is not sufficient to show 
that it can yield S max = 2y/2 for some choice of the parameters. As mentioned earlier, the 
singlet state is completely characterized by the single- and two-particle expectation values. 
Hence, in order to make a comparison with the singlet state, it is necessary to measure or 
compute these two quantities. 



IX. DISCUSSION 



We have presented a computer algorithm that simulates Aspect-type EPRB experiments. 
In the simulation, the source produces particles with opposite but otherwise unpredictable 
polarization. Each particle of a pair is analyzed in an observation station, consisting of a 
polarizer and two detectors (Case I). Placing an additional polarizer in between the source 
and each observation station changes the opposite, unpredictable polarization of the two 
particles into a pair of fixed, but not necessarily opposite, polarizations (Case II). The time- 
tag data of the detection events observed in both stations are used for pair identification. 
Application of quantum theory to both types of experiments yields the single-particle and 
two-particle expectation values that are characteristic for the singlet state (Case I) and the 
product state (Case II). 

The salient features of the simulation model are that: 

• Every essential component of the real laboratory experiment (polarizers, detectors, 
time-tag logic, data analysis procedure) has a counterpart in the algorithm. 
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Identical elements in the experimental setup are represented by identical algorithms. 
For instance, to simulate Case I and II, we use the same algorithm to simulate the 
polarizers. In particular, the algorithm that simulates the polarizer reproduces Malus 
law, which is not essential to reproduce the quantum theoretical results of Case I, see 



the model introduced in Ref. 151 ] . 



• It is event-based and strictly satisfies Einstein's criteria of local causality but it is not 
unique. 

• At any time, it allows free choice of the directions in which the polarization will 
be measured, in contrast to laboratory experiments in which the polarizers in the 
observation stations can take 2x2 directions only 8j. 

• It identifies pairs based on the time-tag of each detection event, using a time window 
W and allowing for several different procedures to define which two photons form a 
pair, just as in real laboratory experiments. 

• To first order in the time window W, it reproduces exactly the single-particle averages 
and two-particle correlations of quantum theory for both Case I and II. 

• It provides information about the frequency of coincidences T. In order to reproduce 
the results of the two-particle correlation as given by quantum theory for Case I, T must 
be sufficiently small. Values of T, corresponding to those found in EPRB laboratory 
experiments [lo| can be reproduced also. For these values of T, the two-particle 
correlation function deviates from the quantum theoretical result but the deviations 
are small. In all cases, the simulation model reproduces the single-particle averages as 
given by quantum theory. 

In our simulation model, the time-tag data are a key element for producing the single- 
particle expectation values and two-particle correlations as given by the quantum theory 
of a system of two 5 = 1/2 spins. In our model, in Case I, the two-particle correlation 
depends on the value of the time window W. By reducing W from infinity to zero, this 
correlation changes from typical Bell-like to singlet-like, without making any change to the 
whole algorithm. Thus, the character of the correlation not only depends on the whole 
experimental setup but also on the way the data analysis is carried out. Hence, from the 
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two-particle correlation itself, one cannot make any definite statement about the character 
of the source. Thus, the correlation is a property of the whole system (which is what 
quantum theory describes), not a property of the source itself. It is of interest to note 
that if we perform a simulation of Case II the single-particle and two-particle correlations 
do not depend on the value of the time window W. In this case, the observation stations 
always receive particles with the same polarization and although the number of coincidences 
decreases with W (and the statistical fluctuations increase), the functional form of the 
correlation does not depend on W. 

Summarizing: We have demonstrated that a simulation model that strictly satisfies Ein- 
stein's criteria of locality can reproduce, event-by-event, the quantum theoretical results for 
EPRB experiments with photons, without using any concept from quantum theory. We 
have given a rigorous proof that this model reproduces the single-particle expectations and 
the two-particle correlation of two S = 1/2 particles in the singlet state and product state. 
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